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Abstract 



> 

C^ ■ To systematically analyze the dynamical implications of the matter content in cosmology, we 

generalize earlier dynamical systems approaches so that perfect fluids with a general barotropic 

^—v ■ equation of state can be treated. We focus on locally rotationally symmetric Bianchi type IX 

^^ \ and Kantowski-Sachs orthogonal perfect fluid models, since such models exhibit a particularly 

^^ ■ rich dynamical structure and also illustrate typical features of more general cases. For these 

fi \ models, we recast Einstein's field equations into a regular system on a compact state space, 

O^' which is the basis for our analysis. We prove that models expand from a singularity and 

t , I recoUapse to a singularity when the perfect fluid satisfies the strong energy condition. When the 

bX)' matter source admits Einstein's static model, we present a comprehensive dynamical description, 

which includes asymptotic behavior, of models in the neighborhood of the Einstein model; these 

results make earlier claims about "homoclinic phenomena and chaos" highly questionable. We 

also discuss aspects of the global asymptotic dynamics, in particular, we give criteria for the 



S . collapse to a singularity, and we describe when models expand forever to a state of infinite 

dilution; possible initial and final states are analyzed. Numerical investigations complement the 
analytical results. 
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1 Introduction 

At present, the matter content in the Universe is a mystery: dark matter, dark energy, quintessence; 
our knowledge that underhes these labels is clearly not as good as we would wish. Apart from ob- 
servational challenges this also poses theoretical issues, for example, very little is known about how 
matter properties influence cosmological evolution, the character of singularities, and the stability 
of various special cosmological models. This motivates the development of a framework that is 
insensitive to the details of the matter content, but still allows one to systematically investigate the 
relationship between properties of matter and cosmological dynamics. The book edited by Wain- 
wright and Ellis P used dynamical systems methods to successfully study perfect fluid models 
with a linear equation of state; in this paper we modify and generalize the formalism in T, and 
attempt to show that the usefulness of dynamical systems methods increases with the complexity 
of the matter source. 

In this paper we consider spatially homogeneous (SH) cosmological models with orthogonal 
perfect fluid matter (i.e., the fluid velocity is orthogonal with respect to the symmetry surfaces) 
described by an effective barotropic equation of state p = p{p), where p denotes the pressure and 
p the energy-density, which we assume to be non-negative; the effective fluid may consist of an 
arbitrary number of non-interacting orthogonal perfect fluids with non-negative energy-densities 
and barotropic equations of state, and a positive cosmological constant. From these assumptions it 
follows that SH Bianchi type I- VIII cosmologies are forever expanding or contracting and compara- 
tively easy to treat, as discussed in the concluding remarks. We therefore focus on SH Bianchi type 
IX and Kantowski-Sachs perfect fluid models which exhibit a particularly rich dynamical structure. 
These models are the only SH models that admit positive intrinsic curvature; they thus naturally 
have a closed spatial topology and we therefore refer to them as closed cosmologies (although all 
Bianchi class A and some class B models also admit closed spatial topologies). However, we re- 
strict ourselves to the locally rotationally symmetric (LRS) orthogonal perfect fluid cases, since 
this suffices to illustrate our approach; generalizations, even to models with no symmetries at all, 
are discussed in the concluding remarks. 

The LRS type IX and Kantowski-Sachs models have been studied before: in [2], 0, |3 dynam- 
ical systems techniques were used to study quite special sources, e.g., perfect fluids with a linear 
equations of state and a cosmological constant; other methods were used in ^5 and pi to treat 
more general matter that included perfect fluids with non-negative pressure, note however that 
this excludes e.g., the inclusion of a cosmological constant. Recently the dynamics associated with 
Einstein's static model was studied in a series of papers that covered several special cases in a LRS 
type IX context: a positive cosmological constant and dust in j7], a positive cosmological constant 
and a perfect fluid with a linear equation of state in [Sj ; the latter work was subsequently general- 
ized to the diagonal non-LRS type IX case in Pj . In these papers it was asserted that the presence 
of Einstein's static model leads to the existence of infinitely many homoclinic orbits [orbits whose 
a- and Lj-limit is the same periodic orbit] that "produce chaotic sets in the whole state space" ; the 
asymptotic configurations of models in a neighborhood of Einstein's model, as of any domain in 
the whole phase space, were claimed to be unpredictable and characterized by "homoclinic chaos." 
We reach a different conclusion: we present a comprehensive description of the dynamics in the 
neighborhood of the Einstein model, moreover, we establish theorems about the global asymptotic 
behavior of solutions. Our results disprove the claims of non-predictability and chaos for models 
close to the Einstein model; furthermore, we prove that the asserted "homoclinic phenomena" must 
be confined to quite small regions of the phase space, if they occur at all; indeed, our numerical 
investigations suggest that "homoclinic chaos" is excluded entirely, but further studies are needed 
to establish this with certainty. 
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Our investigation is based on a reformulation of Einstein's field equations for LRS type IX 
(Kantowski-Sachs) models as a completely regular dynamical system on a compact four-dimensional 
(three-dimensional) state space. Building on this, our analysis of the type IX case rests on three 
cornerstones: methods from dynamical systems theory, in particular center manifold analysis, the 
existence of a conserved "energy" integral, and the close relationship between the dynamics of the 
full system and the Friedmann-Robertson- Walker (FRW) case. 

The outline of the paper is as follows: in Sec. El we reformulate Einstein's field equations for 
the LRS type IX and Kantowski-Sachs models as a regular dynamical system on a compact state 
space; in type IX we also find a conserved "energy." In Sec. El we describe the FRW models; first 
with a potential approach, presumably familiar to a general reader, and then with our dynamical 
systems formulation. In contrast to the potential formulation, the dynamical systems approach can 
be applied to more general cases. In Sec.jllwe investigate the LRS type IX models. Together with 
dynamical systems methods and center manifold reduction we use the conserved quantity to obtain 
comprehensive information about the dynamical possibilities. The dynamical system contains the 
Kantowski-Sachs models as an invariant boundary subset; these models are briefly treated in SecO 
We conclude with a discussion about other cases and how the present work fits in a larger context as 
regards the dynamics of quite general models in Sec. JHl Some proofs and supplementary theorems 
about aspects of the global dynamics are postponed to the appendix, where we also show how our 
variables are related to the Hubble-normalized variables used in e.g., ^. Throughout we use units 
such that the speed of light c and Newton's gravitational constant G are given by c = 1 = 8ttG. 

2 The dynamical systems approach 

The LRS Bianchi type IX models are characterized by the equations 



(1) 



(see p. 40, 123-124 in ^; n2 = n^, o"_ = yields the LRS case). The Hubble scalar H allows one 
to define the length scale i through H = i/i, where i is the cubic root of the volume density. The 
shear is determined by (T+; ni and n2 describe the spatial three-curvature via ^R = —2^1 + 2nin2 
and S+ = n{nin2 — nf). The overdot denotes differentiation w.r.t. clock time t. 

The algebraic equation for p yields H^ + ^nin2 = ^/o + cr^ + j^^i- Since ni and 112 have the 
same sign in Bianchi type IX, we can replace n2 with the variable D := (H"^ + nxn^j'iY''^ ■ We 
then introduce the bounded dimensionless^ variables 

(Qo, Q+) := {H, a+)/D, {M,,nD) := {nJAp) /(UD^) . (2) 

By definition, —1<Qq<1 and Mi > 0, while p > implies 

nD = l-Ql-Mi>0. (3) 

We replace H,a^,ni,n2 with QQ,Q-^,,Ml,D, and p with 0/), which is given in terms of Q^ and 
Ml by (ISJ; instead of p we introduce the dimensionless quantity w = p/p. Then (^ yields one 
decoupled dimensional equation 

D' = -D{QoF + Qo + Q+), F ■.= 2QI - QoQ+ + U'^ + 3w{L))Qd , (4) 



H = 


-H^-2al-l{p + 3p), 


&+ = —3Ha+ — ^5+ , 


hi = 


-{H + 4a+)ni, 


712 = -{H - 2a+)n2 , 


P = 


-3{p + p)H, 


p = 3H^ - 3al + ^^R , 



^For a discussion about dimensions see, e.g., |1U|. 
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and the coupled dimensionless system 



Q'o 


= -(1-Q2)F, 


Q'+ 


= ^M,-l + {Qo-Q+f + QoQ+F, 


M[ 


= 2Mi{QoF-3Q+), 


L' 


= f{L)L{l-L)Qo, 



(5) 



where ' = -^ := jy^ stems from a new dimensionless time variable r; the variable L is discussed 
below. The system ((S| is complemented by the auxiliary equation 

n'D = {2QoF-{l + 3w)Qo + 2Q+)nD. (6) 

The coupled system (jSJ describes the essential dynamics; once solved, the metric can be obtained 
in terms of quadratures, cf. Ch. 10 in ^. The r.h.s. is polynomial in Qo,Q^,Mi. This makes an 
inclusion of the boundary, which consists of invariant subsets, possible; Qq = ±1 yields the forever 
expanding and contracting representations of the LRS Bianchi type II models; Mi = yields the 
Kantowski-Sachs models (see, e.g., [TT) and fi/j = the vacuum boundary. We denote the state 
space of © by 

^X = {(Qo,Q+,Mi,L)}. (7) 

The choice of the variable L depends on the equation of state; L is constructed from i in two 
steps: (i) Let x := i/io, where io = ^(*o) and to is a fixed reference time; introduce X{x) subject to 
the conditions A(x) G ^^[0, oo), A(x) G C^(0, oo), dX/dx > 0, A(0) = 0, and lim2;_^oo A(rc) — > oo. (ii) 
Introduce L := A/(l + A); it follows that L(x) is monotonically increasing, and L £ (0, 1) for x > 0, 
thus L — > describes a singularity while L — > 1 represents a state of infinite dispersion. Moreover, 
we define f{L) := dlnX/dlnx\^(^i^y 

The function u; in Q is to be regarded as a function w{L): from (^ we have dp/{p + p) = 
— 3dx/x, which can be integrated to yield a function p{x) when an equation of state p{p) is given, 
see also the example below. Via x = x{L) this leads to u;(L) := p/p \x(l)- 

The variables ni^2 are algebraically related to the spatial metric and H, o"+ to the extrinsic 
curvature (see Ch. 10 in ^), and are hence continuous in t, or equivalently in x and L; p is 
continuous by virtue of (^, Qo,Q-^-, Mi, D are continuous through 0; however, p, and hence 
w = p/p, is not necessarily continuous - jumps in p correspond to phase transitions. Although 
these can be covered by our formalism, we will for simplicity assume that p{p) is at least C^; 
therefore w{L) and thus the r.h.s. of © is C^(0, 1) (see 12 for related issues). 

For large classes of equations of state, L can be chosen so that © is endowed with the desirable 
differentiability properties for all L £ [0,1], and thus the boundaries L = 0, 1 can be included. 
We define an equation of state to be asymptotically linear if u; — > dp /dp -^ const (= wo,wi) when 
p — > 0, oo. In this case we choose L so that w{L) and f{L) become C^ on L £ [0, 1]; thus w{0) = wq, 
w{l) = wi, /(O) = /o, /(I) = /i. We restrict ourselves to the situation where /o > and /i > is 
possible.^ 

To make things less abstract, let us consider the special case of a source that consists of an 
arbitrary number of non-interacting perfect fiuids with linear equations of state; the total energy- 
density is p{x) = J2i PiO^~^ ^^^^^ where YliPiO = Po (the total energy-momentum T^*^ is a sum of 
the Tj/''^ for each fiuid component i; each component satisfies Tj^'^-i^ = 0, and hence T^'^-i, = 0). We 
assume that pio > and Wi > —1 Vi; causality of each component requires Wi < 1. As an example 
of how to obtain a suitable variable L, let us take the combination of two fluids with wq > wi; 

^By these conditions excessive use of center manifold theory in the dynamical systems analysis is avoided. See |13| 
for a related problem where the general situation is discussed. 



w 
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the choice A(x) = £10 x^w--^!) yields w = wq — {wq — wi)L and f{L) = 3{wo — wi), which clearly 
satisfies the above conditions. As a second example, consider a source that consists of radiation, 
dust and a cosmological constant (a cosmological constant corresponds to pAo = A and wa = — 1); 
a possible choice that satisfies the requirements is A = x, which implies f{L) = 1 and 

(L) = gPrO (1 - L)^ - PAqL^ , . 

^' Pro{l-L)^ + PdoLil-Lf + PAoL^' ^^ 

The dimensionless coupled dynamical system © is thus regular everywhere on the state space 
^X, including the boundaries Qq = ±1, Q^ = ±1, Mi = 0, and O^) = 0; L = 0, L = 1 are also 
included when the equation of state is asymptotically linear, however, even when this is not the 
case, attractors and repellors are distinct sets on L = 0, 1, as we will see in the following sections. 

The system © is invariant under the discrete transformation {Qq,Q+,t) -^ —{Qq,Q^,t), 
which is a consequence of the invariance of the field equations under time reversal. This reflects 
itself in the fact that the fixed points of the system appear in pairs, except for the "Einstein point" 
which is invariant since Qo = Q+ = 0, see below. 

The system (j^l admits the integral 

44/3 n^ mI'^ 

where £' is a constant "energy." This conserved quantity is a generalization of an integral that 
describes the solution structure in the FRW case, cf. Sec.|Sl In ®, the "potential" V is defined as 

V:=-px^, (10) 

where p{x) := p/ po, and hence V <{) since p > 0; y can be viewed as function of x and thus of L 
via x(L), and 

dV/dL = -{l + ?,w)V[f{L)L{l-L)]-'^ and 

V = -(l + 3w)yQo, 

follows from dp/dx = — 3/j(1 + 'w)x~'^ . Eq. (|llj) in combination with Q leads to the claimed 
conservation oi £; £' = 0. The integral £ foliates the interior of the 4-dimensional state space ^X 
into 3-dimensional hypersurfaces. 

As an example of a potential V ^ consider again the combination of two fluids with wq > wi; 
the choice X{x) = {pio/poo)x^'^'"o-'"^'> leads to 

V{L) = -C^ ^-(l+3«,o) (1 _ 2,)l+3«'i \ (12) 

with C = C(poO)Pio) = const. The special case of a cosmological constant and a perfect fluid 
p = wqp results in 

where Vmax = V{^)'-, ^ is defined by w{£) = —3, which in the present case yields £ := {wo + 
l/3)/{l + wo). 

To obtain a better feeling for our formalism we begin with a description of the FRW cosmologies; 
first in the potential approach, presumably familiar to a general reader, and then in the dynamical 
systems picture. 



ViL) 
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3 FRW cosmologies 

3.1 The potential approach 

FRW models are characterized by the scale factor i{t) subject to the three key equations (see, 
e.g., d) 

P = lpf-k, e = -i{p + 3p)i, p = -3{p+p)i/e, (14) 

where t denotes the clock time along the fluid congruence which is orthogonal to the symmetry 
surfaces. The first equation is the Friedmann equation, which is an integral of the second one, 
Raychadhuri's equation; the third equation describes local conservation of energy; k = +1,0,-1 
determines the sign of the spatial curvature. The combinations p + 3p and p + p can be regarded 
as an active gravitational and inertial mass-density, respectively; positive (negative) p + 3p, i.e. 
w > —1/3 {w < —1/3) when p > 0, implies deceleration (acceleration). 

To obtain a dimensionless formulation, i is rescaled w.r.t. a reference time to: x = i/io (with 
io = i{to); see p. 744 of ,14,). A natural choice for a dimensionless time variable is H(to)t, which 
measures time in "Hubble units" ; however, for our purposes we find it more convenient to introduce 
tp := ypo/3 t, i.e., we measure time in "energy-density units." Expressed in x and tp, (fTl)) takes 
the form 

±'^ = E-V , x = i(l + 3w)Vx-'^ , p = -3/5(1 + w) x/x , (15) 

where the overdot now denotes differentiation w.r.t. tp. The equation of state enters via the 
functions w{x) and V{x), cf. (jTU]): note that the potential V is scaled so that ^(1) = —1. The first 
equation in (|15|1 can be interpreted as an "energy" integral 

£ = i^ + y = _^p_3Hyy{pQ£l) = E, (16) 

where the constant E is given hy E = — {3k) / (pqIq) . The problem is thus analogous to the problem 
of a particle with energy E that moves in a potential V{x). The condition V < implies that the 
open (E > 0) and flat models {E = 0) are either forever expanding or contracting; however, the 
qualitative behavior of closed models {E < 0) depends on both the shape of the potential and the 
value of E; we now focus on this case. 

In addition to a non-negative energy-density, we henceforth also assume a non-negative inertial 
mass-density p + p (weak energy condition), i.e., w > —1 when p > (recall that w = —1 
corresponds to a cosmological constant). For simplicity we also assume that there exists at most 
one value of x such that w{x) = —1/3; hence we require dw/dx |^-i(_i/3) < 0, or equivalently 
dw/dp\yj-i(^_il^-^ > 0. Under these assumptions V has a maximum at x = w~^{—l/3), since 
dV/dx = -{l + 3w)Vx-'^ = at X = t(;-^(-l/3) and d^V/dx'^ U-i(_i/3) = -6p^dw/dp\^=_i/s < 0. 

We distinguish five different cases; different matter assumptions yield qualitatively different 
dynamical properties for the associated FRW models, see Tabled the corresponding potentials are 



depicted in Figs. m2(a)l and 2(c 



Figs.^ [2(a)[ and |2(c)| provide a good qualitative understanding of the FRW dynamics, described 
in Table Q however, the exact asymptotic behavior of solutions for x — > 0, co depends on the 
asymptotic properties of the potential. Asymptotically linear equations of state constitute simple 
examples, since w — > wq for x — > and to ^ wi for x — > oo. Consider, e.g., a source that consists of 
non-interacting perfect fluids with linear equations of state; then V{x) = —px"^ = — ^^ pjo2^~ > 

where pio := Pio/po and ^^ pio = 1, and hence limj^^o w = max it; j and lima;^oo w = mintOj. 

3.2 The dynamical system approach 

In the dynamical systems picture the FRW models are represented by solutions of the system © 
on a particular invariant manifold obtained by setting Q+ = and Q', = 0, which implies that 
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Case 


w{x) 


Qualitative FRW Dynamics 


(i) 


w > -1/3, infw > -1/3 


expanding-contracting 


(ii) 


w > —1/3 for small x 

w = —1/3 for a unique x 

w < —1/3 ioT large x 


E = El : forever expanding/contracting 

E = Eg : expanding/contracting-Einstein static 

E = E2 : expanding-contracting or reverse 


(iii) 


w < —1/3, supw < —1/3 


contracting-expanding 


(iv) 


w > -1/3 for X ^ 
and liuix~,ooW = —1/3 


E = El : forever expanding/ contracting 
E = Eg : expanding-Einstein static 
E = E2 ■ expanding-contracting 


(v) 


w < —1/3 for x ^ oo 
and liiUx^ow = —1/3 


E = El : forever expanding/contracting 
E = Es : contracting-Einstein static 
E = E2 ■ contracting-expanding 



Table 1: Matter cases and the qualitative evolutionary behavior for closed FRW models. 



V 




Figure 1: Qualitative shapes of potentials for FRW in the Cases (iii) [small dashed line], (iv) [solid 
line] and (v) [large dashed line] of Table ^ 



Mi = l{l-Ql) (4M( = {l-QlY follows). On the FRW subset the general system reduces to 

Q'o = -U^-Ql){3 + Ql){l + 3w{L)), 



L' 



f{L)L{l-L)Qo 



The conserved energy integral © reads 

£{Qo,L) = 



3(1 -Qg) 
3 + Qg 



V{L) = E<0 



(17) 



(18) 



cf. also (|16() . The conserved quantity determines the orbits of (|17() in the FRW state space {{Qo, L)}. 
To simplify the subsequent discussion, we assume an asymptotically linear equation of state that 
also strictly satisfies asymptotic causality, i.e., wo,wi < 1. 

Case (i), i.e., w > —1/3 everywhere and asymptotically (strong energy condition), cf. Tabled 



yields Fig. 2(a) in the potential approach, while in the dynamical systems approach we obtain 
the flow depicted in Fig. |2(b)[ It is easy to see that the models initially expand, reach a point of 
maximum expansion when Qq = 0, and then recollapse. 
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(a) Case (i) 



(b) Case (i). Phase 
portrait. 





(c) Case (ii) 



(d) Case (ii). Phase 
portrait. 



Figure 2: Qualitative shapes of potentials for FRW Cases (i) and (ii) in Tabled with associated 
flows of the FRW dynamical system ()17() . 



For Case (ii), i.e., matter domination in the limit L ^ 0, (— g < wq < 1) and dark energy 
domination in the limit L ^ 1 (—1 < wi < — |), the potential has the form depicted in Fig. 2(c) 



The associated phase portrait in the dynamical systems approach is given in Fig. 2(d)| the structure 
of the phase space is as follows: 

1. The energy Es = Vmax-, where Vmax is the maximum of the potential, yields the fixed point E 
and the four heteroclinic separatrix orbits, which originate from the fixed points ^F, ^F (or 
E), and end at E (or ^F, \P)- The point E represents Einstein's static universe; its position 
is (0,>C) where C is determined by w{C) = —1/3; ^ F represent flat fluid FRW models with 
linear equation of state (when wi = —1, j_F reduce to the de Sitter model). The separatrix 
orbits divide the phase space into four regions that contain models with distinct qualitative 
behavior. 

2. E = El > Vmax yields regions I and II; region I (II) contains forever contracting (expanding) 
models. 

3. E = E2 < Vmax yields regions III and IV. The orbits in region III (IV) represent models that 
initially expand (contract), reach a point of maximal (minimal) extension at Qq = 0, and 
then contract (expand). 
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The equations in © and H17|l are completely regular, this is in stark contrast to the non- 
regularity of the equations used in [3 |H1 IHl- Non-regularity leads to a misleading picture, e.g., 
compare Fig. 2(d) with Fig.l in "W , where the whole L = subset has been crushed to a point, 
which in turn leads to a deformation of the interior orbit structure; moreover, in contrast to L = 
in our approach, that point is not part of the state space in jH], since the non-regularized equations 
blow up there. 

Of the cases not discussed above, Case (iii) is closely related to Case (i); in Cases (iv) and (v) 
the Einstein point lies on the boundaries L = and L = 1 respectively; although these cases can 
also be treated easily, we refrain from doing so here. 

Thus both the potential and the dynamical systems approaches provide comprehensive infor- 
mation about the qualititative structure of the FRW solution space, however, in contrast to the 
potential approach the dynamical systems approach can be applied to more general problems. 

4 LRS type IX cosmologies 

The foundation of our analysis of the LRS type IX models is the regular dynamical system 0, 
supplemented by the conserved energy ®. For simplicity we assume w < 1 and, initially, an 
asymptotically linear equation of state. We focus on Case (i), where —1/3 < wo,wi < 1, and Case 
(ii), where —1/3 < wq < 1 and — 1 < wi < —1/3. 

We begin with a local dynamical systems analysis: the fixed points and the associated eigen- 
values are given in Table |2 All fixed points correspond to either vacuum solutions or solutions 
that can be interpreted as perfect fluid solutions with a linear equation of state: F represents flat 
FRW solutions for a fluid with a linear equation of state; Q stands for the LRS Kasner solution 
and T for the Taub representation of Minkowski spacetime (see |lj); CS stands for the Bianchi type 
II Collins and Stewart solution; X stands for a solution discussed, e.g., in J^J; finally E represents 
Einstein's static universe. The upper left index refers to L = or L = 1; the lower left index refers 
to Qo = +1 or Qo = — 1, except for ±X where it indicates the signs in the expressions for Qq and 
(5+; the lower right index refers to Q+ = +1 or Q+ = — 1. Note that the fixed points E and ±X do 
not exist in Case (i) while 3^CS do not exist in Case (ii). Table 121 shows that in Case (i) all fixed 
points are saddles except for the source 5}_T_ and the sink ^^.T^. In Case (ii) all fixed points are 
saddles except for '^T_ and LF, which are sources, and ^T-|- and ^F, which are sinks. 



F.P 



Qo 



Qh 



Ml L 



Ai, A2, A3, A4 



^F 
^CS 

±x 

E 



±1 
±1 
±1 

±1 
2 

■(l-3t«i) 







±1 

Tl 

(1+3U.A) 

8 
(l+3wi) 

-(l-Su-i) 











1 

4 



A 
A 
A 
A 
1 
£ 



±(1 + 3wa), ±(1 + 3wa), ^1(1 - wa), ±S{A)fA 
±3{1 - wa), ±2, t4,±S {A) f A 
±3(1 - wa), ±6, ±12, ±S{A)f A 



±i{l + 3wA) 

■ 6(l+3u>i) 



=F 



6(l-t»A) 



=F- 



8 
3 1-wi 



l-3wi ' '2 l-3u)i 



IVb, 



I A iVG. 



{l±i^),±S{A)fA 
iV6 



3uii 



Table 2: Fixed points of the dynamical system © and eigenvalues of the linearizations. Here A 
takes the values or 1; Wa ■= 3(1 - wa){1 + 3wa)/8'^; C is defined by w{C) := -1/3; S{0) = 1, 
S{1) = -1, b := -2w'f{C)C{l -C),c:= {-3w\ + 16wa + 3)/(2(l - wa)) and d := {2Awj + 7wi + 



l)/(l-^i);/(£)>Oand 



-w 



-dw/dL\c > 0, which follows from dX/dx > and dw/dx\c < 0. 
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10 



:q. 



1^ 








(a) —1/3 < w = const < 1 



(b) —l<w = const < —1/3 



Figure 3: Typical phase portraits for —1/3 < w = const < 1 and —l<w = const < —1/3. 

4.1 The w = const case 

In the case of LRS type IX models with a linear equation of state w = const, the equation for L 
decouples from ©; therefore it suffices to consider the state space {{Qq,Q+,Mi)}. We distinguish 
two ranges for w that generate different sets of fixed points: —1/3 < w < 1 and —l<w< —1/3 
(we refrain from giving the special case w = —1/3). Note that the equation system (and thus the 
solution structure) for —l/3<w<l coincides with the system on the subsets L = and L = 1 in 
Case (i) and the subset L = in Case (ii); —l<w< —1/3 is identical to the L = 1 subset in Case 
(ii). The state spaces, and some orbits, are given in Fig. |31 



4.2 LRS type IX models with -| < w < 1 

Next we consider Case (i), i.e., models with —1/3 < w < 1 and —1/3 < wqjWi < 1. We prove that 
all such models expand from a singularity, reach a point of maximum expansion, and subsequently 
recollapse to a singularity. Hereby the cj-limit (a-limit) of every orbit is located on the 2-plane 
(L = 0) n (Qo = -1) {{L = 0) n (Qo = +1)). 

The subcase < t;; < 1 is covered by the general results of jSj and 0, where it is proved, for 
the general Bianchi type IX case, that models do not expand forever, instead they re-collapse to 
a singularity where certain curvature invariants and the density diverge. Our dynamical systems 
approach yields more details about the asymptotic dynamics and the methods of our proof are 
entirely different from those in 0|2j: our starting point is the dynamical system ©. 

On the hypersurface Qo = we obtain F > and therefore Qq < 0; hence, Qo = acts as a 
"semipermeable membrane" for the fiow of the system (jJJ. The hypersurface cuts the interior state 
space in two halves, a future invariant^ half Qq < 0, where L monotonically decreases, i.e., L' < 0, 
and a past invariant half Qo > 0, where L' > 0. It follows that fixed points and periodic orbits are 
excluded in the interior of the state space. 

Application of the principle of monotone functions, see, e.g., 1 , and investigation of the struc- 
ture of the flow on Qo = 0) yields that the w-limit of the future invariant set Qo < must be 
contained in L = 0; analogously the a-limit of every orbit in Qo > must lie on L = 0. More 
specifically, the cj-limit is located on the 2-plane (L = 0) n (Qq = —1). To show this, consider an 
orbit (Qo, Q+, Mi, L){t) in Qo < 0. Since L — > for r ^ oo, the function V, cf. Fig. 2(a) satisfies 



^In Appendix^ we give a brief introduction to dynamical systems nomenclature. 
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V — > — oo. The integral of motion @ fulfills 

(1 - Q2)4/3 ^ D_, ^ Q (r ^ oo) , (19) 

and thus Qq — > — 1 (t — > cxd), as claimed; analogously, the o-limit of Qo > lies on (L = 0) n (Qo = 

As regards the future asymptotic behavior of orbits in Qo > (and the a-limit of Qo < 0), 
there exist, a priori, two possibilities: an orbit either passes through Qo = into Qo < 0, and thus 
converges to (L = 0) n (Qo = ~1)) or it remains in Qo > for all times. However, a proof by 
contradiction shows that the second scenario is excluded: let us assume that the solution satisfies 
Qo > OVr. The monotonicity principle then implies L — > 1 for r — > cxd. Hence, y — > and, 
from (P, 

It follows that either ^d — > or Mi — > as r — > oo. This restricts the w-limit of the orbit to the 
2-surface 

[(L = 1) n (Qo > 0) n {Qd = 0)] u [{l = i) n (Qo > 0) n (Mi = o)] . (21) 

Since the w-limit set of the orbit is located on a 2-surface, the cj-limit is either a fixed point, a 
periodic orbit, a heteroclinic cycle, or it is the whole 2-surface that acts as the attractor. The local 
analysis of +Q+, +T_, and ^F reveals that these fixed points are saddles such that no orbit from 
the interior of the state space can be attracted. Moreover, the surface (j21|) does neither contain 
periodic orbits nor heteroclinic cycles; the possibility that the w-limit of the orbit is the whole 
2-surface is excluded by the structure of the flow as well. Since the 2-surface (|2T|) does not act as 
an (j-limit set, we have arrived at a contradiction; Qo > OVr is thus excluded and therefore the 
w-limit of an orbit in Qo > lies on (L = 0) n (Qo = — 1). 

Note that the "saddle structure" of the surface (PT|) is valid for any value wi > —1/3; it follows 
that it is not necessary to assume asymptotic linearity of the equation of state, i.e., w ^ wi > —1/3 
is not required, (miw) > —1/3 is a sufficient assumption. 



4.3 LRS type IX models with -1 < w < 1 



Consider now Case (ii), i.e., models with —1/3 < if < 1 for L < £, w{C) = —1/3, — 1 < w < —1/3 
for L > C, and w' := dw/dL \c < 0,^ whereby V has a single maximum, cf. Sec. 13.11 

Below we give a comprehensive description of the dynamics of models in a neighborhood of the 
Einstein point and of the asymptotic behavior of these models. The analogous problem has been 
investigated in [THH] for the special case of a fluid with a linear equation of state and a cosmological 
constant, however, we arrive at different conclusions. We demonstrate that the asymptotic behavior 
of models close to the Einstein point is predictable and not governed by "homoclinic phenomena 
and chaos". In addition we also show several global results. 

Our analysis is based on our dynamical systems formulation, i.e., the regular system ((SJ on 
the compact state space ^X; the analysis rests on three cornerstones: the local dynamical systems 
analysis of the Einstein point E, in particular center manifold reduction; the use of the conserved 
quantity £; the recognition and use of the close relationship between the dynamics of the full system 
and the FRW subset. 



For brevity, a prime also denotes difFerentiation w.r.t. the argument L for functions like w{L) and V{L). 
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The first cornerstone in our analysis is the local dynamical systems analysis of the Einstein point 
E. Two eigenvalues, Ai^2 = ±3v6/4, are associated with eigenvectors that lie in the FRW plane, 
vi^2 = (1,0,0, ±2v6/(3ti;')), where b := —2w'f{C)£{l — £); two purely imaginary eigenvalues, 
A3^4 = ibiv6, lead to a rotation in the plane orthogonal to the FRW plane. Hence the linear space 
at E decomposes into the direct sum E^ © E^ (B E'^ oi a one-dimensional stable (unstable) subspace 
£■* (E^), and a two-dimensional center subspace E^ = ((0, 1,0, 0), (0, 0, 1, 0)). 

To adapt to E and the structure E^ © E"^ © E'^ we introduce new variables, 

1/9 w' \ 1-1/9 w' \ 

go = 2 ( Qo + ^ (L - £) J , g+ = Q+ ,mi = Ml - - , I = - IQo - — (L - C) j , (22) 

which transform the dynamical system Q to 



where the nonlinear terms are collected in the respective terms N = N{qQ,q^,mi,l). 

The invariant center manifold M^ is two-dimensional, tangential to -E*^, and contains E. Locally 
it is represented by the graph of a function h : E^ ^> E^ ® E^, 

{q+,mi) ^ {qo,l) = {hg^,{q+,mi),hiiq+,mi)) , (24) 

which satisfies h(0) = and grad /i(0) = 0. Center manifold reduction, see, e.g. |15j . reduces the 
full nonlinear system (|23() to the locally equivalent system 

Al 0\fq,\ fq'A_fO i\fq,\^fN . ^^5) 



I' J yo X2J yij ' ym[j y-^ oj \mij ' yNm^ 

where N = N{{hg^^{q+,mi),q+,'mi,hi{q+,mi)). 

In general it is impossible to obtain the graph of h in terms of elementary functions, however, 
approximate solutions can be obtained by series expansions: 



1/ 9 w'\ , , .rl/ '^ w 



hqoiq+,mi)=qo = -lqo + -Y^ j ^ ^z('?+,"ii) = / = - Uo - g ^ ' (^6) 



where 



go = -2VA;g+mi + 0(||(g+,mi)f), I = ^-hf^kql + prnj + Oi\\iq+,m,)f) , (27) 

where k := -2^/{3w'{2'^ + 36)), b := -2w' f{C)£{l - C). 

The coordinates in ()25() are not the coordinates (|22() . but coordinates adapted to the center 
manifold structure. We assume that the two coordinate systems agree to first order; this is justified 
since E'^ © E^ coincides with the FRW invariant subset and because of ((^ and ((77|) . To avoid 
excessive notation we use the same symbols for the coordinates of (|22|) and (|25j) . 

The reduced system (^^ represents a decoupled system, and hence the properties of the flow in 
the neighborhood of E can be deduced straightforwardly. The center manifold M|; is represented 
by ^0 = ^ = 0, and the induced dynamical system by the (g+, ?7Zi)-system of (HSJ. The solutions are 
"circular" periodic orbits centered about the fixed point E, as we prove below. The fixed point E is 
an (j-limit for two orbits and an a-limit for another pair, which follows from setting g+ = mi = 
and go = (or / = 0); the four orbits are the separatrix orbits in the FRW plane, see Fig. |2(d)| 
Each periodic orbit acts as cj-limit for two one-parameter families of orbits: in the invariant set 
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qo = 0, Eq. ()25() describes one one-parameter family of orbits in the half-space (/ > 0) (another 
one in / < 0) that spiral down (up) towards the periodic orbit at / = 0; similarly, in the invariant 
space I = 0, the periodic orbit is the a-limit for another two one-parameter families of spirals. For 
general initial data {% 7^ 0,/ / 0) the closed orbits act as a "saddles", see Fig. IH] below; we thus 
observe a natural generalization of the FRW picture. 

The second cornerstone in our analysis is the conserved quantity £. To locate E at the origin 
we introduce the variables 

qo = Qo, q+ = Q+, mi = Ml -1/4, l = L-C, (28) 

so that Eq. @ is transformed to 

^^'^°'^+'"'^'^^-(3-4,i-4mi)(l + 4mi)V3-^- ^^9) 

To discuss the 3-hyper surf aces of constant energy, £^ = £', in a neighborhood of E, we expand (|29|) 
up to second order, using V{C + I) = Vmax + V"{C)P where Vmax = V{C) and V" < 0, 

-4— ^go + 2 E +ll(^+'"^i)ll =3|1 ^ ) ' ^^^> 

where ||(g+,mi)|p := 4g+ + 16mf/3. 

We employ £ to prove that the orbits on the center manifold, i.e., the solutions of the (g+, mi)- 
system in (|25j) . are periodic. The idea is simple: the intersection of a hypersurface £ = E with a 
transversal surface that is itself invariant under the flow of the dynamical system yields an orbit 
with energy E on that surface. It is not possible to study M^ D {£ = E) directly, however, the 
investigation of E'^ Ci {£ = E) turns out to be sufficient. We therefore set go = ' = in ((^ . and 
obtain 3 — (3 — 4g^ — 4?ti,i)(1 -|- 4mi)^''^ = 3(1 — Vmax/E), which describes a family of closed curves, 
parametrized by E. In an appropriately small neighborhood of the origin, these curves can be 
approximated by ||((7+,mi)|p = 3(1 — Vmax/E) = const, where E < Vmax, cf. (PU)) with qQ = I = 0. 
Since grad<S ~ (O, cig+, C2(4m,i -|- g^), O) , where ci, C2 are non-vanishing, it follows that the surface 
£" = £^ is orthogonal to E^ at the intersection; this guarantees the intersection to be a closed curve 
also when E^ is continuously deformed to M^. 

Every solution in a sufficiently small neighborhood of E obeys the relation ||((7+,'mi)|p = const 
approximately, since the equations for {q^,mi) decouple and yield ||(g+,mi)|p = const in fl^ . 
Hence the hypersurface ||(g+,mi)|| = const is a (locally) invariant set. Each orbit in this invariant 
set gives rise to a one-parameter family of orbits that differ only w.r.t. the (polar) angle along 
II (g+, mi) II = const. Therefore, one can represent an orbit -modulo its phase- by an orbit in the 
3-space ^X = { (go,^ ||('7+,m,i)|| ) } that fulfills ||(g+, ?tii)|| = const; hereby the four-dimensional 
dynamical picture is reduced to a three-dimensional one. 

The state space ^X is endowed with a Lorentzian scalar product of signature ( — 1-+), Eq. (|3flj) 
describes 2-surfaces of constant Lorentz norm. The light cone, E = Vmax, divides the state space 
into three disconnected regions: E > Vmax yields region I, the chronological past of E, and region 
II, the chronological future; E < Vmax yields region III/IV, the spacelike region, see Fig. [3 Note 
that the disconnected regions III and IV of the FRW picture |2(d)| are merged to one connected 
region III/IV. The coordinates % and / are null coordinates w.r.t. the metric; this can be seen 
from (EU and the relation V"{C) = -3w'Vmax [f{C)C{l-C)-''\-^, which can be derived from (fTT|) . 

The third cornerstone in our analysis is the close connection to the FRW case. In ^X, the FRW 
plane is given by ||((7+,mi)|| = 0. The separatrix orbits, E = 14nax, are given by the null rays 
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to-r\ 



Figure 4: The FRW plane \\{q^ 
Einstein point E as r ^ ±00. 



,mi) 



in the space '^X. Four null rays, 



\i^ 



converge to the 



Qo = and / = 0. Let u^ {u^ ) denote the future (past) null ray along the go-axis, and v^ the null 
rays along the /-axis, see Fig. ^ From Sec. 13.21 we know the global asymptotics, 



(E 
(E 



i-F), 
o_F), 



r,+ 



(°+F - E) , 
(LF^E). 



(31) 



Orbits close to v , are dragged along and thus exhibit similar asymptotic behavior as L -^ 0,1. 
However, in the general picture the FRW attractor '^F is generalized to (L = 0) n {Qq = — 1), which 
describes the contracting Bianchi type II boundary with a linear equation of state characterized by 
wq (the FRW source '^F is generalized to (L = 0)n((5o = +l)i which describes the expanding version 
of the type II boundary); see App. lB.ll A priori the FRW source LF becomes (L = 1) n (Qo = — 1); 
however, on that surface it is only the fixed point LF that can act as an a-limit; a similar a;-limit 
statement holds for "^Y. 

A combination of the ideas of our analysis yields an effective description of the dynamics close 
to E: in ^X, an orbit with a given energy E lies on a hyperboloid (or light cone) of constant Lorentz 
norm, cf. pUjl : simultaneously, every orbit satisfies ||((7+,mi)|| = const. Therefore, the intersections 
of the energy- hyperboloids with the planes ||(g+,mi)|| = const yield a complete description of all 
orbits in ^X; furthermore, the asymptotics of (|31|) and the above generalization determines the 
asymptotics of the orbits. We investigate the scenarios case by case: 




ll(q+. mj)ll 



Figure 5: The light cone E = Vmax in ^X divides the state space in three regions that generalize 



the regions in the FRW picture 2(d) 
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Figure 6: Case E > Vmax- Intersection of the mass hyperboloids of energy E with a plane 
II (g+, mi) II = const in ^X yields two hyperbolas, one in region I and one in region II. The axes are 
the same as in Fig. (jSJ. 



When E > Vmax, Eq. (|3U() describes two "mass" hyperboloids in ^X; one in region I and one 
in region II, see Fig. El Intersection with planes ||(g+,?7ii)|| = const yields a family of spacelike 
hyperbolas which are parallel to null rays asymptotically. Fig. |21 this indicates that the hyperbolas 
of region I originate from LF, as uf does, and end on (L = 0)n(Qo = ~1) li^^s t'". Analogously, the 
hyperbolas in region II originate from (L = 0) Pi {Qq = +1) and end at ^F. The described global 
asymptotic behavior of orbits also remains true for orbits remote from E, as proved in App. El 
Thus the FRW picture is generalized: region I (II) contains models that contract (expand) forever. 

When E = Vmax, Eq. (|3n|) describes the light cone in ^X, see Fig. El Intersection with the FRW 
plane ||(g+, mi)|| = yields the four null rays discussed above, cf. (|31j) and Fig. 01 Intersection with 
planes ||(g+,mi)|| = const > yields a family of hyperbolas that resemble those in the E > Vmax 
case; the global behavior is thus analogous. 

Finally, when E < Vmax, Eq. ()3U() describes one hyperboloid of positive Lorentz norm in the 
region III/IV in ^X. We distinguish three subcases: 

Case (a): Intersection with planes ||(g+,?TT.i)|| = const < 3(1 — '^°'' ) yields a family of timelike 
hyperbolas, see Fig. 8(a) This leads to a generalization of the behavior of orbits in region III and 
IV in FRW: we distinguish models that originate from LF (resp. (L = 0) n {Qq = +1)) and end in 

^F((L = 0)n(Qo = -l)). 

Case (b): Intersection with a plane ||(g+,?TT.i)|| = 3(1 — 'g'"' ) yields the point ||(g+,mi)|| = 
3(1 — '^"^ ) , Qq = 0, I = 0, which represents a periodic orbit, and four in/outgoing null rays 
(cf. also the discussion following Eq. H3U() 1. see Fig. |8(b)1 Two of these null rays represent orbits 
that originate from either (L = 0) n (Qo = +1) or LF and converge to the periodic orbit as r — > oo; 
the other two originate from the center manifold orbit and end at either (L = 0) n {Qq = —1) or 
^F. The FRW picture is thus generalized; periodic center manifold orbits (one for each value of 
E) take the place of the Einstein point E, see Fig. [71 This figure also indicates how the abstract 
^X picture is related to the '^X state space picture, cf. Fig. 8(a 

Case (c): Intersection with planes ||(g'+,?TT.i)|| = const > 3(1- 
like hyperbolas, see Fig. 8(c)[ this case resembles the E > Vmax case and yields similar behavior. 



and Fig. 
'"""' '' results in a family of space- 



Accordingly, there exist models that expand forever from a singularity to infinite dispersion, despite 
that their energy E is lower than the "potential barrier" Vmax- In this case ||(g+,?7ii)|| acts as an 
additional energy that is sufficiently large to lift the orbit over the top of the hill from the "side of 



recollapse" to the "side of expansion" in Fig. 2(c) when it passes through the neighborhood of E. 
When E is sufficiently close to Vmax , and if the orbit passes through a sufficiently small neighbor- 
hood of E, the presented arguments about global asymptotic behavior are mathematically rigorous. 
The known asymptotics of the null ray orbits (|31|) , and continuous dependence on initial data, guar- 



0,li 



antees that orbits that are close to the null rays are dragged to neighborhoods of j! F. The fixed 
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(a) Spiral tubes that de- 
scribe the stable manifold 
of a periodic orbit; on each 
side of the periodic orbit 25 
spiral orbits are depicted. 



(b) Spiral tubes that de- 
scribe the unstable man- 
ifold of a periodic orbit; 
each side of the periodic 
orbit consists of 25 spiral 
orbits. 



(c) The four spiral tubes 
of Figs, (a) and (b) in the 
state space; the diameter of 
the tubes is too small to be 
resolved. Note that in the 
projection of Figs, (a) and 
(b) the tubes would have 
overlapped. 



Figure 7: Illustration of the close relationship between the dynamics close to E and FRW dynamics. 
At the scale of Fig. (c) there is no visible difference between this figure and the analogous figure 
depicting the FRW separatrices in the same projection, except that the orbits do not converge to 
^F in general, but to the generic limits ^T-t, cf. Add. IbTTI and Fig. IIUI In this figure the matter 
content is radiation and a cosmological constant, however, the qualitative behavior is the same for 
any matter source that admits the Einstein static solution. 



point ^F is a sink, hence all orbits that pass through a sufficiently small neighborhood of \F end in 
^F. An analogous statement holds for the source LF. When an orbit is dragged along D~ towards 
'^F, it represents a model that contracts when L < £; in App. Elwe show that it must continue 
to contract and reach (L = 0) n (Qo = ~1)- Analogously, orbits that head for "^F along Uj^ (for 
T -^ — oo) expand forever from (L = 0) n (Qo = !)• 

We emphasize that the dynamics of the LRS type IX models is closely related to the FRW 
case and completely predictable near the Einstein point E as regards asymptotic states: in the 
FRW case, the asymptotic states are determined by the energy E; in the LRS type IX case, the 
asymptotic states are determined by the energy E and by the value of ||(g+,?TT.i)|| = const. This 
disproves the claims of j3 Ej about non-predicability of asymptotic states for models close to E. 
Moreover, as shown in Appendix^ the homoclinic behavior of orbits asserted in [7||H], can be 
disproved not only for orbits that pass through a neighborhood of E, but also for large regions 
of the state space, such as L < C, L > Ce- in these regions the asserted orbits whose a- and 
w-limit is the same periodic orbit cannot exist. Furthermore, numerical results reveal that the 
qualitative features described above for a sufficiently small neighborhood of E carry over to a large 
neighborhood of E and a large range of energies, as can be seen in Fig. Inland Fig. 1101 This indicates 
that the assertions about "homoclinic phenomena and chaos" are wrong for the entire state space. 

Finally note that the discussion of the dynamics in the neighborhood of the Einstein point E is 
independent of the assumption of asymptotic linearity of the equation of state and that the global 
asymptotic behavior of orbits remains valid in a wide sense. The boundaries are no longer part of 
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Figure 8: Case E < Vmax- Intersection of the timelike hyperboloids with energy E with different 
planes IKg+jTTii)!! = const. We obtain timelike hyperbolas in Case (a), a vertex and a null cone in 
Case (b), and spacelike hyperbolas in (c). The vertex of the null cone in (b) represents a closed 
orbit in the center manifold. The global asymptotic behavior of the orbits is governed by the null 
ray orbits i>j , cf. Fig. 0] The axes in the figures are the same as in Fig. Q. 

the systems, however, it is still correct to state that the a- and w-limits of orbits are located on 
(L = 0) n (Qo = ±1) and (L = 1) n (Qo = ±1)- 

5 Kantowski- Sachs cosmologies 

The coupled system of equations for the Kantowski-Sachs models is obtained from the general 
system © by setting Mi = 0; we thus have the three-dimensional state space {{Qo,Q+, L)} 
depicted in Fig. 1111 The coupled system admits a number of equilibrium points whose associated 
eigenvalues are given in Table 01 



F.P 



Ai 



A. 



A. 



^F 



± 



X 



±{1 + 3wa) 
±2 
±6 

I l-3«)i 



^1(1 -wa) 
±3(1 - wa) 
±3(1 - wa) 



-1-3 1-wi 
'2 1-3W] 



-1-3 1-wi 
'2 l-3u;i 



±SiA)fA 
±SiA)fA 
±SiA)fA 

1-Vd 



Table 3: Fixed points and eigenvalues for the Kantowski-Sachs submanifold. Here A = 0, 1, 
S{0) = 1, S{1) = -l,d= {2Awl + 7wi + 1)/(1 - wi). 



It follows that in Case (i) all fixed points are saddles except for ^T_, which is a source, and 
'^T^., which is a sink. In Case (ii) all fixed points are saddles except for ^T_ and LF, which are 
sources, and *'_T+ and ^F, which are sinks. 

In Case (i) all models expand from a singularity, reach a point of maximum expansion, and 
then recollapse to a singularity, as follows from an analysis similar to that of the LRS type IX case. 
The situation in Case (ii) is more complicated; local analysis limits the possibilities somewhat but 
numerics is required to give a detailed picture. Fig. 1111 illustrates various possibilities, e.g., there 
are solutions that expand from singularities to infinitely dispersed isotropic states and solutions 
that contract from infinitely dispersed isotropic states to singularities. 
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(a) Periodic orbits act 
as "saddle" orbits in the 
neighborhood of E. 



(b) The center manifold 
dominates even far away 
from E. 



(c) The same orbit as in (b), 
now depicted in the state 
space through a {Qo,Q+,Mi)- 
projection. 



Figure 9: Example of the effect of the center manifold on dynamics. The source is orthogonal dust 
and a positive cosmological constant, however, the pictures depict behavior that is typical for all 
sources that admit a unique Einstein point E. 



6 Conclusion and outlook 

In this paper we have regularized Einstein's field equations and obtained a dynamical system on a 
compact state space for the LRS Bianchi type IX and Kantowski-Sachs perfect fluid models with a 
barotropic equation of state. This allowed us to systematically study the effects of matter properties 
on cosmological evolution, and we obtained a global picture of the dynamical possibilities. In type 
IX, the cornerstones of our analysis were: center manifold theory in connection with Einstein's 
static solution -when existent; an "energy integral" that severely restricted the dynamics; a close 
connection between the FRW and LRS type IX cases. This led to a comprehensive description of the 
dynamics in the neighborhood of Einstein's static model: given an initial value in the state space in 
the neighborhood of Einstein's model, it is possible to predict whether the solution will expand to an 
infinite dilute final state or eventually contract to a singularity. This disproves previous assertions 
that the dynamics near the Einstein model exhibits non-predictability and chaos. The present 
formalism is generalizable to cover more general models such as general Bianchi type IX models for 
which it is also possible to derive a conserved energy which, together with center manifold theory, 
again severely restrict the dynamics; indeed, we expect that the conclusions drawn in this paper 
will have direct bearing on this problem as well. In addition to this we have obtained several global 
results; again we expect our methods and results to be generalizable to the general Bianchi type 
IX case. 

The present approach is also of relevance for the forever expanding Bianchi type I- VIII perfect 
fluid models. As in the present case one can introduce a variable L and obtain an autonomous 
formulation, however, in this case one can also use ln(£/^o) ^-s a time variable (the natural time 
variable when H is used as the normalization factor, as discussed in [J), which then leads to 
that w becomes a time dependent function and thus that one obtains a non-autonomous problem. 
However, in the case of an asymptotic linear equation of state one obtains an asymptotically 
autonomous problem and for certain problems one can then apply the Strauss- Yorke theorem to 
determine the asymptotic dynamics ^Hl, ^Zj- Both approaches are expected to have advantages 
and disadvantages that depend on the problem at hand. Although one can draw the conclusion 
that Bianchi type I- VIII perfect fluid models with p > are forever expanding, the details of how 
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(a) Spiral tubes that emanate from 
a periodic orbit far from the Ein- 
stein point. On each side of the 
periodic orbit, 25 orbits with the 
same energy that describe the sta- 
ble manifold of the periodic orbit, 
are depicted. 



(b) Spiral tubes that emanate 
from a periodic orbit far from 
the Einstein point. On each side 
of the periodic orbit, 25 orbits 
with the same energy that de- 
scribe the unstable manifold of 
the periodic orbit, are depicted. 



(c) The four spiral 
tubes of Figs, (a) and 
(b) in the state space; 
they emanate from the 
stable Priedmann/Taub 
points on the expanding 
and contracting Bianchi 
type II boundary 
subsets. 



Figure 10: Illustration of the dynamics far from the Einstein point: the qualitative structure of 
the flow is still the one of a small neighborhood of E. Clearly there is no "homoclinic chaos." In 
this figure the matter content is radiation and a cosmological constant, however, the qualitative 
behavior is the same for any matter source that admits the Einstein static solution. 

this actually happens is a quite intricate problem due to asymptotic self-similarity breaking, which 
occur for the more general classes of models |18| , ^7j ; we expect asymptotic self-similarity breaking 
to hold for the same class of models also for asymptotically linear equations of state, but the details 
will depend on the equation of state. 

We thus expect that some of the ideas in this paper are generalizable to other SH models, 
however, we also believe that some ideas are generalizable to models with fewer symmetries, even 
to models with no symmetries at all, so-called Go models, particularly when it comes to asymptotic 
dynamics associated with singularity formation. Again one can treat a barotropic equation of state 
by introducing a function similar to L that parametrize the equation of state suitably. In special 
cases of matter dominated singularities additional assumptions like asymptotic linearity will have to 
be imposed on the equation of state, however, generically the singularity is expected to be "vacuum 
dominated" ^01; in this case we expect that one can treat equations of state with quite general 
asymptotic behavior. 
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A Dynamical systems nomenclature and definitions 

Consider a dynamical system y = f{y) and the associated flow ^f The w-limit uj{x) (a-limit a(x)) 
of a point x is defined as the set of all accumulation points of the future orbit {$t(x) 1 1 > 0} (past 
orbit {^t{x) 1 1 < 0}) of x. The w-limit of a set S is uj{S) = \J^^gUj{x). An orbit {^t{x) | i G E} is 
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Figure 11: The Kantowski-Sachs submanifold. Case (ii). 



called homoclinic, if there exists a (fixed) point z such that a{x) = uj{x) = {z}; an orbit is called 
heteroclinic if it originates from a fixed point zi and ends in a different fixed point Z2- A set S is 
called future (past) invariant, if $t(5) C S'Vt > (Vt < 0). 

The monotonicity principle states that if there is a function Z on an invariant set S that is 
strictly decreasing along orbits, then uj{S) = {s £ S\S\ liuix-^s Z{x) ^ sup^Z}, and analogously 
for the Q-limit. 

Two dynamical systems are called (locally) equivalent, if there exists a homeomorphism (of 
coordinates) ^ such that the flows are topologically equivalent $j = ^ 
conditions the equivalence is a diffeomorphism. 



o $^ o $. Under certain 



B Asymptotic states and global dynamics 

In this section we discuss aspects of the global dynamics of case (ii) of the LRS type IX models: first 
we characterize the asymptotic states of models that converge to a singularity or to an infinitely 
diluted state, then we describe the global qualitative behavior of models with E > Vmax, and 
thirdly we treat some aspects of the global dynamics of models with E < Vmax] the considerations 
complement the investigation of Sec. 14.31 which focussed on the dynamics near the Einstein point 
E. 



B.l Asymptotic states 

To investigate the asymptotic behavior of orbits that converge to a singularity as r — > oo we make 
use of the conserved energy Q. In analogy with (fT9|) we find Qo — *■ — 1 as r — > oo, i.e., the attractor 
is the set (L = 0)n (Qo = -!)• In the FRW picture |2(d)1 the attracting surface (L = 0)n (Qo = -1) 
reduces to the fixed point ^F; in the general picture (L = 0) n {Qq = — 1) contains four fixed points: 
'^F, 'LQ_, '^T-i-, and '^CS. The generic attractor is ^T+, a three-parameter family of orbits ends 
there; ^F attracts a two-parameter set of orbits, '^CS a one-parameter set, cf. Table |2j Since no 
periodic orbits or heteroclinic cycles exist on (L = 0) n (Qo = — 1) (see jTj), the fixed points are the 
only possible w-limits. Analogously, the surface (L = 0) n (Qo = 1) contains the a-limit of orbits 
that expands from a singularity; the treatment of the fixed points is analogous. 

When the equation of state is not asymptotically linear, i.e., when w -/^ wi, but infw > —1/3 
as L ^ 0, then L = is not included in the state space and the local analysis of the fixed points 
breaks down; however, since ^T+ is a sink for all —1/3 < t^i < 1, the point (—1, +1,0,0) remains 
a generic attractor of a three-parameter family of orbits also in the general case. 
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The asymptotic behavior of orbits that expand to a state of infinite dispersion is quahtatively 
simpler. As r — > oo the orbit approaches the attracting set (L = 1) n (Qo = 1); to see that we use 
the same arguments as above. However, on (L = 1) n (Qo = 1)) among the fixed points ^F, +Q+, 
5,_T_, it is only the point ^F that attracts orbits from the interior of the state space, see Table |21 
Analogously, LF is the only source on the set (L = 1) n (Qo = ~1)- These statements remain valid 
also when the equation of state is not asymptotically linear. 

B.2 E > Vmax 

We prove that models that satisfy E > Vmax, i-e., models in regions I and II, are forever expanding 
or forever contracting between a singularity and a state of infinite dilution. Consider a solution 
{Qo,Q+,Mi,L){t) with E > Vmax and consider Q divided by Vmax, 



^ _ 3 ^-1., ^2x4/3, . -1/3 ^(-^) 



^-n\l-Ql)'^'M['/';^. (32) 



V 44/3 D V N^oy 1 Y 

^max ^ ' ^max 

Since V {L) /V^ax ^ 1 and by assumption E/Vmax < 1 it follows that 

^,il-Qlf"<^DMl'^<±,, (33) 

where O/jM-j^ has been bounded by its maximal value. We therefore conclude that Qq > 
const^ Vr, hence either Qo(r) > const or Qo(t) < —const Vr. From the system (jSJ it follows 
that L' > const f{L)L{l — L) or L' < —const f{L)L{l — L) and the claim is established. 

Let us juxtapose the above argument with the discussion in Sec. 14.31 Recall that E > Vmax 
corresponds to a pair of mass hyperboloids in the neighborhood of the Einstein point E. Thus, 
locally, it follows that Qo (= Qo) is either always positive or always negative when E > Vmax, hence 
L' < or L' > 0; the above argument proves that the mass hyperboloids extend to hypersurfaces 
that satisfy Qo > or Qo < globally in the state space. 

Models that satisfy E = Vmax are forever expanding or forever contracting between a singularity 
and a state of infinite dilution; however, five models are exceptions: the Einstein model and the 
four FRW models that correspond to the separatrix orbits - the latter are forever expanding or 
contracting but converge to E. To prove the claim consider solutions (Qo, Q+, Mi, L)(r) with 



E = Vmax', from (|n^ we obtain 



K 



^(l-Qg)^/^ = l^z.My^:^<^ (34) 

for all L y^ C. {^dMI^^ attains its maximal value 3/4^/^ at (Q+,Mi) = (0,1/4).) Eq. (jSU 
implies Qo(''") ^ and hence L'{t) ^ when L{t) / C. Assume that there exists tq such 
that L'{tq) = 0. It follows that L(tq) = C and Qo(''"o) = 0; via the equation in 1)34(1 we obtain 
(Q+(ro),Mi(ro)) = (0,1/4), i.e., (Qo, Q+,Mi,L)(to) = E, hence the solution coincides with the 
Einstein model. Therefore, except for the Einstein universe, all models with E = Vmax are forever 
expanding or contracting (so that L{t) is strictly monotonic). Assume that a model satisfies 
L{t) ^ Lo 7^ 0, 1 for r ^ ibcx); then L'{t) -^ (at least for a subsequence (T„)„gN) and Qo(''") -^ 
as T ^ ±00; from ()3l|) we derive that Lq = C and Q+{t) -^ 0, Mi(t) -^ 1/4. Therefore, the model 
has E as an a- or w-limit, so that it must correspond to one of the four FRW separatrix orbits, 
cf. Sec. 14.31 We conclude that apart from the five special cases all models with E = Vmax satisfy 
L — > 0, 1 for T — > ±00, and the claim is established. 
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B.3 E < Vmax 

Models with E < V^ax exhibit more intricate behavior as aheady seen in the discussion of Sec. I4.3t 
however, chaotic dynamics is excluded from the outset in large parts of the state space, as we will 
see in the following. 

Models that contract for some L < C continue to contract until they reach a singularity. 
Similarly, models that expand for some L < C must have expanded from a singularity. The proof 
is analogous to that in Sec. 14.21 when L < £, the surface Qo = is a semipermeable membrane, 
Q'o < on Qo = 0- Thus, (Qo < 0) n (L < £) is a future invariant set with V < 0, while 
(Qq > 0) n (L < £) is a past invariant set with L' > 0; taking the structure of the flow on Qo = 
into account, the claim follows. When L = C, Q'q < still holds; in fact one can show that Q'q < 
when Q+ / 0; when Q+ = we have Q'q = Qq = 0, but Qq = —4(1 — 4Mi)^, which is negative 
everywhere except for at the Einstein point. Hence Qq = is semipermeable for L < C and the 
claim is established. 

From the above it follows that all models that pass through Qq = when L < C are expanding- 
contracting, i.e., these models expand from a singularity, reach a point of maximum expansion, and 
re-contract to a singularity. 

A model with energy E < Vmax that expands for some L > Ce continues to expand forever to 
a infinitely diluted state. Similarly, a model that contracts for some L > Ce must have contracted 
throughout its past history. Here, Ce is uniquely defined by V{Ce) = E, Ce > C; note that 
Ce — > C when E -^ V^max- The proof of the claim is a variation of the arguments used in IB. 21 
consider the hyperplane Qq = 0; in Qq = 0, the surface £ = E is given by the set of all (Q+, Mi, L) 
such that E/V{L) = (3/4^/3) f^-^Mj'^''^ For each L = const with C < L < Ce this equation 
defines a closed curve centered about (Q+, Mi) = (0, 1/4); for increasing L -^ Ce, the curves shrink 
to the central point (0, 1/4). For L = const > Ce the intersection {Qq = 0)ri{£ = E)ri{L = const) 
is empty. Hence, an orbit with energy E cannot cross Qo = when L > Ce, which implies that 
(Qo > 0) n (L > Ce) is a future invariant set with L' > and (Qq < 0) n (L > Ce) a past invariant 
set with L' < 0. From this the claim follows. The statement can be strenghtened considerably: 
there exists Ce < Ce such that the statement remains valid when Ce is replaced by Ce- The 
quantity Ce is defined as the minimum of all L = const such that F < (and hence Qq < 0) 
on (Qo = 0) n (<f = £') n (L = const) (for F see Eq. Q)); in Fig. El L is chosen according to 
Ce > L = const > Ce'- hence (Qo = 0) n (£" = £^) n (L = const) is not empty, but it is contained in 
the region F < 0. By definition, for all orbits with energy E, for all L > Ce, the hyperplane Qo = 
acts as a semipermeable membrane: orbits can only pass through Qq = from Qo < to Qo > 0. 
It follows that (Qo > 0) n (L > Ce) is a future invariant set with L' > and (Qo < 0) n (L > Ce) 
a past invariant set with L' < 0, and hence the claim is shown. A very good approximation for the 
value of Ce is given by the solution of the equation V{Ce) = E (4/3)(l — w^Ce))"^- 




Figure 12: In Qo = 0, for L = const > Ce, the curve £ = E is entirely contained in the region 
F < where Q'q > 0. 
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C H-normalized variables 

We here give a comparison between the variables used in the paper and the important H-normahzed 
quantities used in e.g., ^ (note that in contrast to our variables ii^- normalized variables break down 
when H = 0): 

Ql ' 4MiQ2 ^ ^ 

^^^_^__^_^__^ g ^ 2S^ + 1(1 + 3w)nH = ^ ^^2°^+ . (36) 
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